Biophysical Chemistry

Intrinsically disordered proteins (IDP) or proteins with intrinsically disordered regions (IDR) perform a plethora of functions mostly in a cellular environment. As unfolded proteins, IDPs can adopt molten globule or coil ensembles of conformations. Regarding the latter the question arises whether they are describable as a self-avoiding random coil. Locally, this requires that amino acid residues sample the entire sterically allowed region of the Ramachandran plot with very similar probabilities and independent on the conformational dynamics of their neighbours. However, various lines of experimental and bioinformatic evidence suggest a more restricted, side chain and nearest neighbor dependent conformational space for individual residues. Over the last 25 years short peptides and coil libraries were employed to determine conformational propensities of amino acid residues in unfolded states. The question arises whether conformational ensembles obtained from these two sources are comparable. In this paper, a variety of metrics were used to compare Ramachandran plots of a limited number of GXYG peptides (X,Y: guest residues) with XY dimers in the coil library of Ting et al.(PLOS 6, e1000763, 2010). The results reveal major differences between corresponding plots, which might in part due to the fact that solely the influence of one of the two neighbours of a given residue is probed by the above coil library while averages were taken over the respective opposite neighbours. The presented results suggest that coil libraries alone might not be a sufficient tool for determining the characteristics of statistical coils of IDPS and IDRs alike.

Mutations in MAPT, the tau gene, give rise to forms of frontotemporal dementia and parkinsonism linked to chromosome 17 (FTDP-17T), with abundant filamentous tau inclusions in brain cells. Some mutations that encode missense and deletion variants can give rise to a clinical picture of Picks disease and filaments made of three-repeat tau. Here we report the electron cryo-microscopy (cryo-EM) structures of tau filaments from individuals with MAPT mutations D252V, G272V, S320F and {Delta}G389-I392. The two-layered Pick fold was present in the individuals with mutations D252V and {Delta}G389-I392. By contrast, mutations G272V and S320F gave rise to a more open variant of the Pick fold, with residues 272-341 rotated by 20-25{degrees} with respect to the rest of the structure. These findings show that missense mutations within the filament core can modify the Pick fold, generating closely related structural variants. In addition, we were able to reconstitute the Pick fold and some of its variants using seeded assembly with recombinant 0N3R tau carrying 12 serine or threonine to aspartate substitutions (PAD12) and missense mutations D252V, G272V or S320F. This work provides a foundation for the development of structure-based diagnostic and therapeutic approaches.

The KRAS oncogene, central to cellular signaling via MAPK and PI3K-AKT pathways, is a notorious cancer driver frequently activated in pancreatic, colorectal, and lung carcinomas. Regulation of human KRAS oncogene expression is important due to its capital role in cell growth, proliferation, and survival. Misregulation of its expression contributes directly to the development and progression of multiple types of cancer. In previous studies, the role of G-quadruplexes elements in both the promoter and 5 UTR regions have shown to play important roles in KRAS expression, particularly when these G4s elements interact with regulatory protein hnRNPA1. In this study, we reveal that KRAS expression is also modulated at the post-transcriptional level through the formation of RNA G-quadruplexes (rG4s) situated at the 5 untranslated region (5UTR) of the mRNA. Biophysical and binding studies were carried out to probe the interaction. Through isothermal titration calorimetry (ITC), we quantified a strong binding affinity between the UP1 domain of hnRNPA1 and short-nucleotide RNA segments capable of adopting different G-quadruplex fold. The binding interaction is characterized by a favorable Gibbs free energy change in the range of {Delta}G {approx} -32 to -34 kJ/mol, suggesting a specific and energetically favorable association. One-dimensional and two-dimensional 1H-15N HSQC NMR spectroscopy revealed pronounced chemical shift changes in residues of both RNA recognition motifs (RRMs) of UP1, signifying direct contact with the rG4 structure.

The importance of human aldehyde oxidase (hAOX1) has increased over the last decades due to its involvement in drug metabolism. Inhibition studies concerning hAOX1 are extensive and a common reducing agent, dithiothreitol (DTT), was recently found to inactivate the enzyme. However, in previous crystallographic studies of hAOX1, DTT was found to be essential for crystallization. To surpass this concern another reducing agent used in crystallization trials. Using tris(2-carboxyethyl)phosphine (TCEP), a sulphur-free reducing agent, it was possible to obtain well-ordered crystals from hAOX1 wild type and variant, hAOX1_6A, which diffracted beyond 2.3 [A]. Instead of the typical star-shaped crystals of hAOX1, at pH 4.7, plates are obtained in the orthorhombic space group (P22121) with two molecules in the asymmetric unit. Activity assays with the enzyme incubated with both reducing agents show that contrary to DTT, TCEP does not lead to irreversible inactivation of the enzyme. The replacement of DTT with TCEP in crystallization of hAOX1 provides a strategy to circumvent enzyme inactivation during crystallographic studies, allowing future applications of new assays, such as time-resolved crystallography.

Proteins with intrinsically disordered regions (IDRs) migrate at a higher apparent molecular weight in sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) complicating their analysis and identification. Here, we investigate the sequence determinants of the hypomobility of IDRs using a series of synthetic low complexity domains. We find that negative charge increases the apparent molecular weight, but neutral polar tracts also have abnormally slow migration. Positive charge and hydrophobic residues decrease the apparent molecular weight, although lysine residues show a biphasic effect with decreased migration at high fractional contents. Combinations of residues show that different sequence contributions to the apparent molecular weight are not additive. The results can be rationalized by the protein-decorated micelle model by considering both SDS binding and the compaction of protein SDS-complexes.

Magic-angle spinning (MAS) solid-state NMR (SSNMR) has been widely used to determine amyloid fibril structures at atomic resolution. Such studies typically rely on homogeneous fibril preparations that produce narrow linewidths and high spectral resolution, enabling reliable resonance assignment and structural analysis. However, many biologically relevant amyloid aggregates are structurally heterogeneous, resulting in spectral broadening and reduced sensitivity that hinder atomic-resolution characterization. Lipids are known to modulate amyloid aggregation pathways and promote the formation of toxic species that are often less homogeneous, further complicating NMR-based investigations. Here, we evaluate the feasibility of utilizing the benefits associated with high-field (1.1 GHz) SSNMR for studying ganglioside GD3-catalyzed A{beta}42 aggregates. Uniformly-13C,15N-labeled A{beta}42 was incubated with GD3 to generate lipid-associated aggregates and analyzed under MAS conditions. 13C cross-polarization magic-angle spinning (CPMAS) spectra and 2D 13C-13C chemical shift correlation experiments using CORD (COmbined R2nv-Driven) mixing were acquired and compared with data collected at 600 MHz. Despite the heterogeneous nature of the GM1-associated assemblies, the 1.1 GHz spectra exhibit enhanced sensitivity and improved spectral resolution. Better resolved resonances corresponding to selectively structured regions of A{beta}42 are observed, indicating the presence of an ordered core within the lipid-associated aggregates. These results demonstrate that ultrahigh-field SSNMR significantly improves the characterization of heterogeneous amyloid assemblies and provides a promising approach for atomic-level investigation of biologically relevant, lipid-modulated A{beta} aggregates.

Multi-domain proteins (MDPs) adopt diverse conformations arising from cooperative inter-domain motions, and such dynamics are intimately coupled to their biological functions. Quantitative characterization of these motions is crucial for elucidating their functional mechanisms. Although small-angle X-ray scattering (SAXS) provides information on overall domain arrangement, the limited experimental constraints hinder reliable discrimination of conformational ensembles derived from molecular dynamics (MD) simulations. To address this limitation, complementary experimental constraints that enable to observe domain-selective structural information are required. Inverse contrast-matching small-angle neutron scattering (iCM-SANS), combined with segmental deuteration, enables selective visualization of individual domains and thus provides such complementary information. However, practical strategies for preparing segmentally deuterated MDPs with arbitrary domain labelling have yet to be established. Here, we develop an experimental protocol that integrates controlled protein deuteration with high-efficiency multi-step protein ligation to generate a segmentally deuterated MDP in high yield. The combined use of SAXS and iCM-SANS yields complementary structural constraints that enhance discrimination of MD-derived conformational ensembles. This protocol expands the applicability of segment-selective visualization and also provides an opportunity for high-precision analysis of dynamics in complex MDPs. SynopsisSegmental deuteration enabled by high-efficiency multi-step protein ligation, combined with inverse contrast-matching SANS and SAXS, provides structural constraints that improve discrimination of molecular dynamics ensembles of multi-domain proteins. IMPORTANTthis document contains embedded data - to preserve data integrity, please ensure where possible that the IUCr Word tools (available from http://journals.iucr.org/services/docxtemplate/) are installed when editing this document.

Protein-DNA complexes are involved in vital cellular functions like gene regulation, replication, transcription, packaging, rearrangement, and damage repair. In this work, streamlined geometric formalism for computing the absolute binding free energy was used to obtain chemical accurate in silico estimation of binding free energy of three Protein-DNA complexes. Additionally, molecular interactions between Protein and DNA involved hydrogen bonds, electrostatic, van der Waals, and hydrophobic interactions. Using this formalism, researcher can obtain the absolute binding free energy for a Protein-DNA complex with remarkable accuracy and modest computational cost.

Deoxyribonucleic acid (DNA) stands as one of the most foundational concepts in life sciences, essential for students to master. However, when surveyed about the forces that stabilize the double-stranded DNA structure, many students exhibited a conceptual bias-- favoring base pairing as the primary stabilizing force, while overlooking the equally critical role of base stacking interactions. To investigate the origins of this misconception, students conducted a comprehensive analysis of 35 widely used textbooks. Their findings revealed that one-third of these texts explicitly emphasized base pairing as the sole stabilizing force in their written content. Furthermore, two-thirds of the textbook contained illustrations that reinforced this bias, visually highlighting base pairing while neglecting base stacking. Recognizing this bias, students embarked on a literature review to gain a more accurate and nuanced understanding of DNA stabilization. Through this research, we identified three concept areas--DNA structure and function, environmental effects on DNA, and DNA-protein interactions--to illustrate how base pairing and base stacking work in concert to stabilize the antiparallel double helical structure of DNA. This interplay between base pairing and base stacking is crucial not only for the structural integrity of DNA, but also for its biological functionality. By addressing this conceptual bias, we aim to promote a more balanced and scientifically accurate representation of DNA stabilization in educational materials.

The structural and dynamic properties of membranes are known to vary with bilayer size and hydration. While Molecular Dynamic (MD) simulations are a powerful tool for studying cellular membrane systems, the results can be sensitive to the analysis work-flow and software. In this study, all-atom MD simulations (500 ns) were conducted on systems of 256, 512, and 1024 POPC lipids at 40, 80, and 160 waters per lipid. With these simulations, a two-fold study was performed: (1) to assess the convergence of structural and dynamic properties of POPC bilayers as a function of membrane size and hydration level using CPPTRAJ (CPP), including area per lipid (APL), bilayer thickness, order parameter, headgroup orientation, and lateral diffusion, and (2) to compare the analysis output and performance of four software packages: CPP, GROMACS (GRO), MDAnalysis (MDA), and LiPyphilic (LiP). For the first objective, our results show that the average values of the bilayer thickness, order parameter, and headgroup orientation are largely independent of the size and hydration levels studied. In contrast, lateral diffusion coefficient was sensitive to both size and hydration. We found that increasing the system size primarily decreased the statistical variance of the APL and thickness. For the second objective, all four packages produced consistent results for APL and thickness, with the most significant discrepancy being a known artifact from the gmx order tool when applied to unsaturated carbons. Performance bench-marks identified CPP as the fastest serial tool for all properties, whereas parallelization benefited MDA and LiP in some metrics. These findings provide a practical roadmap, demonstrating that moderately sized systems (e.g., 256L), combined with an optimized tool such as CPP, offer an efficient workflow for membrane structural property analysis.

Granulocyte macrophage-colony stimulating factor (GM-CSF) is a cytokine that plays a role in immune modulation. Its expression is associated with a multitude of different effects ranging from harmful, as in diseases such as rheumatoid arthritis and multiple sclerosis, to beneficial, as in the case of mitigation of diabetes type I and neutropenia. However, there is a large gap in knowledge explaining how GM-CSF toggles its structure for such physiological and pathological interactions. Our work describes an ongoing attempt to address this gap by focusing on a clustered histidine triad within -helices near the N-terminus, which prior studies have suggested play a role in binding ligands at an acidic pH. While GM-CSF is known to be highly flexible at a more acidic pH, several properties of its histidine triad remain unclear at the physiological pH at which GM-CSF would encounter its binding partners. We describe an effort to characterize the role of the GM-CSF histidines under physiological pH, specifically to determine if these histidines are key to GM-CSF structural integrity, and whether individual histidine residues modulate binding as they do at a lower pH. Our findings reveal that, while the histidine residues have an impact on GM-CSF structure, flexibility, and stability, they alone do not modulate the affinity for ligands at neutral pH. These data provide an initial explanation for the pleiotropic functions and interactions of GM-CSF within a biophysical context. Graphical Abstract O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=82 SRC="FIGDIR/small/700583v1_ufig1.gif" ALT="Figure 1"> View larger version (24K): org.highwire.dtl.DTLVardef@a6fffcorg.highwire.dtl.DTLVardef@1f00c30org.highwire.dtl.DTLVardef@b04c50org.highwire.dtl.DTLVardef@6224d9_HPS_FORMAT_FIGEXP M_FIG C_FIG

Insulin amyloid aggregation is a key pathological and pharmaceutical concern, particularly in the context of Type-2 Diabetes (T2D), where amyloid deposition of protein can impair therapeutic efficacy and contribute to cell death leading to local tissue damage. Although gangliosides--glycosphingolipids containing sialic acid residues--are known to modulate amyloid formation in neurodegenerative disorders, their influence on insulin aggregation remains largely unexplored. In this study, we investigate the effects of gangliosides GM3 and GD3 on insulin aggregation. Using Thioflavin-T (ThT) based fluorescence kinetics, Fourier Transform Infrared (FTIR) spectroscopy, Circular Dichroism (CD) spectroscopy, Small Angle X-ray Scattering (SAXS), Nuclear Magnetic Resonance (NMR) spectroscopy, and Transmission Electron Microscopy (TEM), the aggregation pathway, changes in the secondary structure and morphology of insulin aggregates have been characterized. Our results show that both GM3 and GD3 lipids accelerated insulin aggregation in a concentration-dependent manner while steering the pathway away from classical fibril formation, producing short, beaded structures distinct from the extended fibrils observed under lipid-free conditions. CD and FTIR data analyses revealed that insulin in the presence of gangliosides formed non-fibrillar intermediates with distinct secondary structures: {beta}-sheet-rich globular clusters in presence of GD3 and -helical intermediates in GM3-treated samples. Cytotoxicity assays further demonstrated that ganglioside-induced aggregates are significantly less toxic to cells when compared to insulin-only aggregates. Furthermore, ganglioside-bound insulin oligomers retain seeding capacity, suggesting that they can nucleate further aggregation despite their non-fibrillar morphology. These findings underscore the role of gangliosides in modulating insulin amyloid polymorphism and toxicity, offering new insights into their potential impact on the pathology of T2D and treatment strategies. O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=121 SRC="FIGDIR/small/703542v1_ufig1.gif" ALT="Figure 1"> View larger version (28K): org.highwire.dtl.DTLVardef@5bf40eorg.highwire.dtl.DTLVardef@f400ddorg.highwire.dtl.DTLVardef@164dcd8org.highwire.dtl.DTLVardef@def4e7_HPS_FORMAT_FIGEXP M_FIG C_FIG HighlightsO_LIGangliosides GD3 and GM3 accelerate insulin aggregation, forming non-fibrillar assemblies. C_LIO_LIGanglioside-bound insulin aggregates are less cytotoxic than fibrillar aggregates. C_LIO_LIDespite altered morphology, ganglioside-bound aggregates retain seeding ability. C_LI

Amyloid beta (A{beta}), one of the hallmark proteins of Alzheimers Disease (AD), aggregates into plaques that are strongly linked to cognitive decline and neuronal death. Reducing its aggregation propensity may provide a strategy to slow the progression of AD. While chirality modulation has emerged as an innovative approach to disrupt this process, research has primarily focused on alterations at the C position, often overlooking the impact of the second chiral center, such as the C{beta} atom of Threonine. Furthermore, the underlying mechanisms governing these chiral effects remain elusive. Given the intrinsically disordered nature of the A{beta} peptide, we employed temperature-replica exchange molecular dynamics (T-REMD) simulations to explore its rugged conformational landscape. We considered sequence mutations (A2T, A2V), N-terminal chirality inversion of the first six residues (A2V1-6D and WT1-6D), and alteration of the second chiral center (C{beta}) of Threonine (A2TC{beta}). By analyzing the effect size and population change induced by these mutations and chiral modulation, we concluded that the modulation at the N-termini is not confined locally but also exerts specific effects on the central hydrophobic core (CHC) region. Inspection of their free energy landscape and representative structures reveals that the protective or pathogenic effects of these variants correlate with their similarity to the wild type (WT) ensemble. Beyond these static thermodynamics analyses, a direct connection to phase transitions was made by estimating heat capacity as a function of temperature. Both analyses predict that A2TC{beta} may exert a pathogenic effect, in contrast to the protective nature of A2T. These findings offer a deeper understanding of the effects of site-specific mutations and chirality and shed light on the development of advanced therapeutic strategies for AD.

Alzheimers disease (AD) is a leading cause of death among the elderly, with no existing treatment. The development of therapy is further challenged by a limited understanding of molecular pathogenesis and the absence of reliable early-detection biomarkers. Neuroimaging and lipidomic studies reveal structural and biochemical alterations in both gray and white matter in AD patients, including disruptions in membrane organization and neuronal signaling pathways. In the present work, we employed lipidomics-guided modeling of membranes in gray and white matter regions in healthy and diseased (AD) conditions, and used all-atom molecular dynamics (MD) simulations to examine how AD-associated alterations in lipid composition influence the structure, spatial organization, and micro-heterogeneity of neuronal plasma membranes in different brain regions. Data suggest that Alzheimers disease-associated lipid alterations in gray matter (GM) and white matter (WM) impact membrane thickness and microdomain distribution, highlighting the critical role of lipid composition in maintaining neuronal membrane homeostasis and function. Higher-order cholesterol-ceramide- sphingomyelin-enriched domains are more abundant in the neuronal membranes of the GM region in diseased conditions. Under AD-mimicking conditions, lipidomic analyses demonstrate that neuronal membranes in GM experience more substantial compositional and structural remodeling than those in WM. Our results show significant changes in membrane microdomain distribution across the lipid bilayers, and, interestingly, these changes are more pronounced in the gray matter than in the white matter. This study establishes a framework for modeling the tissue-specific lipidomics data to understand how disease-driven compositional changes affect the structure, organization, and dynamics of biological membranes.

Filament-forming proteins such as TasA (Bacillus subtilis) and camelysins CalY1, CalY2 (Bacillus cereus) pose a particular challenge for structural analysis due to their strong tendency to self-association and their polydispersity, which severely limits their ability to crystallize or to be a target for NMR-spectroscopy. To address this, it is necessary to modify the amino acid sequence to prevent filamentation. Engineering a series of N- and C-terminal truncated variants by removing flexible parts is often key to success. N-terminal extensions are also a powerful tool for obtaining crystals of fiber-forming proteins.

We have conducted all atom molecular dynamics simulations of POPC and DPPC lipid bilayers using AMBER Lipid21 force field with eight different water models, including SPC/E, TIP3P, TIP3P-FB, TIP4P-FB, TIP4P-Ew, TIP4P/2005, TIP4P-D, and OPC, to identify the most compatible one without any modification. A number of parameters have been computed in order to understand the structure of the lipid bilayer: Area per lipid, Isothermal compressibility modulus, average Volume per lipid, electron density profile, bilayer thickness, X-ray and neutron scattering form factors, deuterium order parameter, and radial distribution function. The estimated Area per lipid, Isothermal compressibility factor, volume per lipid and bilayer thickness are highly consistent with experimental results for the SPC/E water model, indicating its suitability with the AMBER Lipid21 force field, insted of any modification. The bilayer electron density profiles of both the lipid bilayers demonstrate a little augmentation of water penetration with respect to the membrane surface for TIP4P-D water model. However, the experimental X-ray and neutron scattering form factors are aligning well with the simulated results for all studied water models, and TIP4P-D shows better for X-ray data. The deuterium order parameter for lipid acyl chains value less than 0.25 for all observed water models, depicting their disorderness for both the lipid bilayers. The lateral diffusion and reorientation autocorrelation function of the lipid molecules in both the bilayers are computed to reveal their dynamics across all water models. In comparison to other water models, the simulated trajectories predict better structure and reasonably fair dynamic properties for the SPC/E water model. The TIP4P-Ew water model reproduces the lateral diffusion co-efficient in close agreement with experiment. Reorientational dynamics for both the lipids in the bilayers for eight different water models are observed; the presence of slow and slowest time components corresponds to the lipid axial motion (wobble motion) and Twist/Splay motions. So, in view of the overall performance of the different water models with the AMBER Lipid21 all atom force field in reproducing membrane physical properties, the SPC/E water model appears to be an optimal choice.

Protein stability arises from a fine balance between stabilizing forces such as hydrophobic interactions, hydrogen bonding, and ionic interactions, and destabilizing contributions from solvent exposure and electrostatics. Although hydrophobic burial is the dominant driving force for folding, intra-chain hydrogen bonds and ionic interactions modulate stability in context-dependent ways, with effects that vary depending on their location and environment within the protein. Most studies of protein stability have focused on perturbations induced by pH, solvent composition, or mutations in protonated water, leaving the influence of solvent isotopes relatively underexplored. Notably, despite stronger hydrogen bonding in D2O, proteins exhibit diverse stability responses upon transfer from H2O to D2O, suggesting that differential hydration of nonpolar groups plays a key role. Here, the solvent isotope effect on protein stability is investigated using double-chain monellin (dcMN), a {beta}-sheet-rich, two-chain protein with well-characterized folding behavior. By combining conventional equilibrium unfolding measurements with hydrogen-deuterium exchange mass spectrometry (HDX-MS), the stability of wild-type and a less hydrophobic mutant (C42A) dcMN was compared in H2O and D2O, revealing greater stabilization of the wild-type protein in D2O and highlighting the importance of hydrophobic interactions in governing isotope-dependent stability.

Autophagy is a critical neuroprotective mechanism, the impairment of which can lead to severe neurodegenerative diseases. Spinocerebellar ataxia type 1 (SCA1) is a monogenic neurodegenerative disorder, characterised by the presence of protein aggregates and consequent loss of cellular functions. The expression of mutant Ataxin1 (ATXN1) in glial cells has been demonstrated to induce inflammatory responses and loss of supportive functions, thereby exacerbating neuronal degeneration in SCA1. Autophagic dysfunction has been shown to affect both neurons and glial cells, resulting in widespread pathological consequences. In this work, we aimed to evaluate the efficacy of two small-molecule autophagy activators, AUTEN-67 and AUTEN-99, in models of glia-specific SCA1 in Drosophila. Our results demonstrate that AUTEN-99 has a stronger autophagy enhancing effect, with significantly improved response times and survival rates, compared to untreated ATXN1 mutants. Glia-specific assays in mouse primary hippocampal cultures also confirmed that AUTEN-99 is a more effective activator. Ultimately, co-treatment of neuronal and glial cultures did not reveal any synergistic benefits from combining the two AUTEN compounds compared to single-agent treatment. Our findings contribute to a better understanding of the utility of AUTENs and may help to understand the critical role of autophagy in neurodegenerative diseases.

Intercalation of small molecules between DNA base pairs affects DNA conformation, disrupting essential cellular processes including replication, transcription, and repair. We investigated conformational changes in 18-mer DNA upon intercalation of doxorubicin, SYBR Gold and YOYO-1 using extensive MD simulations. Two main patterns for the intercalation were identified: RISE-type intercalation occurs between adjacent base pairs and extends the DNA helix with decreased twist angles, while OPEN-type intercalation proceeds through base-pair opening without significant DNA extension. Kinetic analysis revealed that association rates for intercalation followed the order: first YO-moiety (mono-intercalation) > SYBR Gold > doxorubicin > YOYO-1 (bis-intercalation). Free energy landscape showed that forces at DNA termini reached up to 117 pN during stretching. Notably, base pairs adjacent to intercalators were protected from strand separation, accompanied by additional helical unwinding. MM-PBSA/GBSA analysis revealed that the driving force for intercalation is the stacking energy, and the binding affinity was highest for minor groove binding. Persistence length decreased with single molecule binding but recovered with two molecules due to their electrostatic repulsion. Mechanical properties of intercalated DNA showed position-dependence, demonstrating that multiple intercalation modes coexist in solution. The heterogeneous nature of intercalation explains why experimental measurements reflect ensemble averages rather than single binding configurations.

Ligand binding has been shown to induce significant alterations in the conformational landscape of proteins. Traditional crystallography approaches have provided valuable input about the end states in ligand-binding reactions. However, dynamical relationships between ligand binding and backbone rearrangement often remain obscured by crystallographic structures. In the present study, we use time-resolved serial synchrotron crystallography (TR-SSX) to directly visualize indole binding in the cavity of T4 lysozyme L99A in microcrystals under controlled environmental conditions. By integrating fixed target crystallography with LAMA-based ligand delivery, we have been able to track the progression of ligand binding and backbone rearrangement. By utilizing an occupancy refinement protocol, we have been able to quantify structural populations. Our studies reveal that ligand binding for this protein cavity follows a diffusion-limited process that progressively rearranges the F -helix of the protein towards a dominant conformational state. These findings establish an observable link between ligand diffusion, occupancy evolution and conformational adaptation within a crystalline environment. More broadly, our work shows how TR-SSX can quantify ligand and conformational populations during binding, providing a framework to interpret structural adaptation in real time.